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.-nSTRACT 

An  attempt  has  been  made  to  analyze  several 
factors  v;hich  contr-ibute  to  the  difference  of  motion 
betv/oen  an  ascendir.f-;  balloon  and  acti.il  v/ind  profiles, 
Tv.-o  analytical  solut-ions  are  presented  v/hich  give  the 
v.'ind  response  error  of  the  balloon  as  a  function  of 
the  altitude,  the  v;ind  profile,  and  the  location  of 
o/.e  balloon  \\fithin  the  wind  shear  layer.   These 
solutions,  both  including  2nd  neglecting  the  apparent 
nassj  agree  viell   v.'it,h  computer  solutions.   The  apparent 
mass  of  ac;:ual  balloons  has  been  measured. 
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I.   INTRODUCTION 

The  satisfactory  recording  of  magnitude  and  direc- 
tion of  wind  at  Vc rious  altitudes  by  means  of  an  ascending 
balloon  depends  among  other  things  upon  the  capability  of  the 
balloon  to  assume  the  velocity  of  the  surrounding  air  in  the 
shortest  possible  time.   The  difference  between  the  actual 
air  movement  end  the  respective  motion  of  the  ascending 
balloon  is  called  "response  error"  (Ref  1),  and  a  balloon 
which  would  assume  instantaneously  the  actual  air  velocity, 
would  have  a  "zero  response  error." 

Physically  the  phenomenon  of  response  error  can  be 
described  as  follows,  '.ihen   the  ascending  balloon  enters  into 
a  layer  of  horizontal  wind  which  Jte  s  a  nearly  constant  velocity 
versus  altitude,  the  balloon  has  to  be  accelerated  until  the 
relative  velocity  between  bf;lloon  and  surrounding  air  becomes 
insignificant.   During  the  period  of  acceleration,  we  will 
have  a  decreasing  response  error  which  a.pproaches  zero  when 
the  balloon  velocity  is  proCtically  equ^tX- to  the  wind  velocity. 

If  an  ascending  balloon  enters  into  horizontal 
wind  which  increases  in  magnitude  with  altitude,  the  balloon 
has  to  be  steadily  accelerated  while  it  passes  through  the 
wind  layer.   This  results  in  a  response  error  depending  on 
the  variction  of  the  wind  v;ith  altitude. 

These  simple  consider.:  t  ions  indicate  that  a:  balloon 


with  a  high  response  error  may  lead  to  faulty  information 
about  the  wind  pcttern,  and  in  the  following, an  attempt 
will  be  made  to  analyze  sevrrcl  reasons  which  contribute 
to  the  response  error. 

The  movement  of  the  ascending  balloon  is  governed 
by  the  equation  of  motion  which  represents  essentially 
Newton's  Second  Law.   A  closer  study  of  this  equation 
indicates  that  a  balloon  will  follow  the  actual  air  move- 
ments rapidly  if  its  speed  of  ascent  is  low,  its  ratio  of 
weight  to  surface  area  is  small, and  if  its  apparent  mass 
is  a  minimum. 

Because  of  certain  conveniences,  spherical  balloons 
are  mostly  used  for  this  purpose  and  this  selection  establishes 
certain  aspects  in  regard  to  drag  coefficient,  Reynolds 
nu-ber  effect,  and  the  magnitude  of  the  apparent  mass. 
This,  however,  does  not  mean  that  spherical  balloons  are 
really  the  best  objects  for  this  purpose. 

The  question  of  how  .light  a  balloon  per  unit  of 
cross  section  area  can  be  built  is  a  problem  of  balloon 
fabrication.   The  a pparent  mass  of  a  spherical  balloon  is 
proportional  to  the  third  power  of  the  radius  and  since  the 
apparent  mass  tends  to  increase  the  response  error,  the 
smallest  balloon   ossible  should  be  used.   However,  the  bal- 
loon has  to  be  large  enou -h  to  climb  to  the  required  altitude 
in  the  prescribed  time  and  must  also  be  visible  on  the  radar 


screen.   Therefore,  the  entire  matter  of  designing  the  most 
suitable  bclloon  is  on  optimization  process  in  which  the 
weight  per  unit  of  projected  areb ,    the  apparent  mass,  tlie 
drag  coefficient,  snd  the  balloon  buoyancy  are  the  principal 
components. 

The  pursuit  of  this  total  problem  exceeds  the 
scope  of  this  study,  which  shall  be  limited  to  an  analysis 
of  the  trajectory  and  to  a  measurement  of  the  a pparent  mass 
of  actual  balloons. 

The  trajectory  analysis  shov/n  in  Section  II  is 
pursued  in  a  closed-form  solution  and  a  first  approximation 
both  with  and  without  the  effect  of  appsrent  mass.   The 
results  of  these  efforts  are  compered  with  earlier  publica- 
tions (Ref  2).   An  additional  calculation  is  then  performed 
which  includes  the  term  of  the  apparent  mass  of  a  sphere 
derived  from  potential  flov/  theory. 

In  the  third  section, measurements  of  the  apparent 
mass  of  actu?.l  balloons  are  described.   The  influence  of  the 
Reynolds  number  upon  the  apparent  mass  appears  to  be  important 
since  the  magnitude  of  the  apparent  mass  is  a  function  of 
the  flow  pattern  related  to  a  subcritical  or  supercritical 
Reynolds  number  operation.   The  import'ence  of  this  investi- 
gation is  emphasized  through  the  fact  that  most  ascending 
spheres  used  for  the  detection  of  wind  operate  in  the  super- 
critical, transitional,  and  subcritical  Reynolds  number  regions 


-II.   TR/iJECTORY  CALCULATIuIJ 

A.       Development  of  the  Equations  of  Kotlon 

The  manner  in  which  an  ascending  balloon 
follows  the  actual  horizontal  movement  of  air  .and- vlind 
follows  from  the  equation  of  motion  of  the  balloon, 
Newton's  second  law  may  be  v^-ritten  symbolically  as: 


tF  --  m,,,^  .  (1) 


where 


m^  .  ■=  total  mass  of  the  balloon 
tot 

V    =  absolute  velocity  of  the  balloon 

IF    ^   sum  of  all  external  forces  acting 
on  the  balloon. 

V/e  choose  Z  as  the  vertical  and  X  as  the 
horizontal  direction.   The  motion  of  the  balloon  is 
assumed  to  be  restricted  to  the  X-Z  plane. 

The  total  mass  (m^^^)  of  Eqn  1  consists  of: 

1  )  mo  =  mass  of  the  balloon  material,  including 
valves,  hardware,  etc, 

2)  Tcu     =  mass  of  the  gas  within  the  balloon 

3)  m'   =  apparent  mass. 

The  apparent  mass  arises  from  the  transfer  of  kin- 
etic energy  from  the  shpere  to  the  surrounding  air  and  for 
ideal  fluid  has  been  found  to  be  equal  to  half  the  mass  of  ai; 


displaced   by   the   balloon    (Ref  3). 

The   extern.  1   forces   acting  on   the   balloon  are 
aerodynamic   drag,    gravity,    end    buoyancy   forces.      Mathemati- 
cally,   this    may   be   expressed   as: 

IF-D+W^B^  '(2) 

where 

D  =  aerodynamic  drag 

17  =  weight  of  balloon  material  plus  included  gas 

B  =  buoyancy. 

We   must    here   differentiate   between  three  different 
velocities.      The  absolute   velocity  of  air   is  denoted   by   V^, 
that   of   the   balloon   by   V,    and  the    relative   velocity   of  the 
balloon  v/itli   respect   to  the   air  by  V    .      By  this   definition 
of   V     we    have : 

^~-  V- v;,  •  (3j 

Writing  the  velocities  in  terms  of  components,  we  have 

^      A      A 

V  -  V^  I  +  V^  k 

A  (4) 

We    have    here   assumed   that   the   wind   lias   no   vertical   component. 
Thus,    we   find    from   Eqns   3   and    4   that: 

Vr   ^(V,-V^)i   .v/k    .  (5) 

V/e   are   nov;  in  a    position  to    v/rite  expressions  for 
the    external    forces  acting   on    t!ie   ualloon.      The   drag  force 


is  conventionally  expressed  as: 


A 


D^-Cd^PVp'sT  (6) 


where 


•D2 


Cp  =  coefficient  of  drag 


p   ■=  air  density 

S  =  cross  sectional  area  of  balloon  (also  the 
projected  area) 

A 

^     =  unit  vector  in  the  direction  of  the  relative 
velocity  (V^), 

A 

and    X    may  thus   be   expressed  as: 

consequently   we    find: 

D  -CoipVrSV;  .  (8) 

The  gravity  force  is  given  by: 

W  =-(  m3.nn^)g1l  .  (9) 

Finally,    the    buoyant   force    is   given  by  the   weight   of   air 

displaced   by  the   balloon.      Thus: 

_  A 

B  =pVo|  Kg  ,  (10) 


where 


Vol   =  volume    of  balloon 

g        =   acceleration   of    •-'■ravity   =   32.2   ft/sec^ 

k        =  unit   vector   in  the   Z  direction. 


V/e   therefore   obtoin  the   resultant   of  the   forces: 

2f.-Cd^P  SV,  \^  -(m3^nQ.)gk^p  gVolk    .  ^^^^ 

Using   relation    (5)    in    (11)    we   find: 
F  =  C^pSV,[(V,-Vji-V,k]-(nn3.m^)gk-  pgVol^k   .  (12) 


Substituting  this   rel^ition,    along  with  Eqn  k,    into   Eqn  1 
yields: 

-(nn^+m^)gk  +  p  Vol  g'k 
Or,    in  terras  of  the   component   directions; 


(  m,  .  m,^.  m' )( V,^  -  V,  k  )  =  -Coip^  V[(  V^--  Vj  i  -  \4  k] 

A        ,    A 


(13) 


(m^+m^^m')V;^=  -  C^^p  S'^i'^-^yJ  '        (14) 

(  m^-m^+m')  V^^  -CpipSVpV2-(m^^rT^)g+ p  g  Vol  ,       (15) 

where  the  magnitude  of  Vn  is  given  by:  ■ 

V.  ^(^^-Mcv^'^V/'  .  (16) 

These  equations  correspond  to  Eqns  9a  end  9b  of  Ref  2,  the 
only  difference  being  that  Ref  2  did  not  consider  the  effects 
of  apparent  mass.   Since  the  apparent  mass  constitutes  a 
fairly  large  percentage  oi'  the  physical  mass,  it  is  ne.cessary 
th6t  this  be  taken  into  account  for  accurate  results.   A 
comparison  of  the  solution  which  includes  apparent  mass  with 
one  where  it  is  neglected  is  presented  in  Section  C. 


The  mass  of  gas  within  the  balloon  can  be  ex- 
pressed as: 

m^^  PgVoI  ,  (17) 

where  p^    °  density  of  the  gas  within  the  balloon. 

We  now  v;ill  proceed  to  represent  the  internal  gas 
density  {  Pq)    in  terms  of  the  external  air  density.   The 
internal  gas  density  is  given  by  the  perfect  gas  law  as: 

P^  -  R,T,  •  .  (18) 

V/e  assume  that  a  pressure  sensing  valve  is  used 

such  th^t  tl'ie  difference  between  the  internal  and  external 

pressures  is  kept  constant,  or: 

p-  p  rAP  =  constant  ■  - 

or 

P^  r  P  ^AP  •  (19) 

Substituting  Eqn  19  into  Eqn  l8  yields: 

Pr~-    -F^    ^^^      ■  (20) 

Assuming  that  the  balloon  is  in  thermal  equili- 
brium with  the  surrounding  air,  we  may  assume  that  TrT=T^4Lj,=T 


and  write: 


p  =  -^  .  AE  .  .  (21) 


g 


Since   we   nay  write    (for  the   ambient  air): 
P^pRT     , 


we   obtain: 


p  J_R_p^[  AP 


o  G 


(22) 


(23) 


Substituting  Eqn  23    into  Eqn  1?  gives: 


m 


^R^ 


G 


\%l 


pVol.  + 


AP 


TR, 


Vol 


(24) 


The  apparent  mass  is  conventionally  represented  as: 

nn'=KpVol  (25) 

(Note:   for  potential  flow,  we  have  K  -=  |.) 

Utilizing  relations  (24)  and  (25)  in  (14)  and  (15) 
yields: 

(m^-  -R-Pvoi  .^VoUKpVoDV^  --   S£^(v_-V,  )     (Ua) 


s   R 


TR 


XW   X 


(m,.  ^PVol.  ^P-Vol^KpVol,)V  -  -  ^^?r^^  V  - 


R 


TR 


2 


-(m^.^-p  vouA£^Vol)g.pgvol  . 

These  may  be  rewritten  in  terms  of  design  parameters  as; 


(15a) 


2 


XW   X 


(26) 


(Q^qp  )  V,  =-  --^-^ -  (Ci.g-Volp)g.pgvol,  (27' 

G 


2 


where: 


Ci 


-s  ^  ff,  vol 

R 


C2--  (-g-  *K  )   vol     . 

We    shall   attempt   to   obtain  a    closed-form  solution 
for  these    equations    in  the   following   section. 


B. 


Solutions  of  the  Governinp;  Equations 


We  shall  examine  two  solutions  in  this  section. 
The  first  shall  be  a  more  or  less  "exact"  solution  wliich 
accounts  for  the  variation  in  air  density  over  the  altitude 
interval  under  consideration.   The  second  will  assume  a 
constant  air  density  over  this  interval  and  proves  to  be 
more  easily  evaluated.   Finally,  the  two  solutions  will  be 
compared  and  some  general  results  shown. 

We  will  here  concern  ourselves  with  shear  layers 
which  are  not  extremely  thick,  e.g.,  altitude  intervals  less 
than  $00  ft.   For  this  size  of  interval,  we  may  safely 
assume  that  the  vertical  velocity  remains  constant.   Thus, 


V. 


0  and   v;e   may   solve   Eqn   2?   for   V     and   obtain 


R 


Vr 


P   (1 


R 


AP 


Vol  g  -  m^  g 


(29) 
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Substituting  this  relatioh  into  Eqn  26  yields: 


(C.1.C2P)  v^  =  ^-^^^^^— ^(C3P -Ci  )g    '  (30) 

where : 

C3^(l-|-JVol     .  (31) 

V/e   may   expand   V,,-   as: 

v    -    dVx       dVxdp  _   dVx  dP    dj2     ^  ,      . 

x"dt'dpdt"dpdhdt 

One    recognizes  that  ~=  V^..      Thus,    we   have: 

0     -vi^dVx       .  (33) 

^x    -  ""^  d  h  d  p 

d  P 
In  order   to   evalxiatc   — —     ,    v/e   make   use   of   the    exponential 

d  h 

density  lav/  given  by:  '  ''  ' 

£  =  e'  (34) 


Po 


> 


where: 


P  =  air  density  at  altitude  h 

Po  =  air  density  at  altitude  h   (h^  is  the 
altitude  of  beginning  of  shear  layer) 

P  =■  logarithmic  density  slope. 

Differentiating,  we  obtain: 

dp  -/3(h-ho) 

^=-Po/3e^     -/3p  .  (35) 

Substitutinc;;  this  relation  in  Eqn  33  we  find: 

Vx  '-PP"^^'-     ■  (36) 
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utilizing  this   equation,    relation    (30)    becomes,    after 
rearranging,: 


lYx      (  H;w-  Vx) 
dp    '     f3pV| 


Ci-qp 


(37) 


In  view  of  the   objective   to  find   the  difference 
betv/een  the  wind  and   balloon  velocity    (wind    response    error)  ^ 


it    is   convenient   to    introduce 


d  Vxw 
dp 


at  this  time.   Subtract- 


ing this  quantity  from  both  sides  of  Eqn  37  yields: 


d(M,-Vxw)  _  (M<w-Vx) 


dp 

or,  rearranging: 

d(\^v-^)-  g 


PPV; 


crC3P 


_    d  Vxw 
y        dp 


dp        ■  PPvf 

This  eaui-tion  has  the  form; 


C1-C3P 


(V-VJ 
xw    >^ 


dV; 


xw 


p 


(38) 


dV 
dx 


f^x)  y  =  g(>^) 


(39) 


The    integrating   factor   for  the   left    side    of   tliis   equation 
is   given  by: 

e  .  (40) 

Multiplying   both   sides   of  Eqn  39  by  Eqn  4O   gives: 
rf(x)dx    ^v  /f(x)dx  ;f(x)dx 

^  7n7  -  e  f(x)y  =  g(x)e  .  (41) 


e 


d  X 


One    recognizes   tliat 


6_ 
dx 


ff(x)dx     l  /f(x)dx       y 


rr(x)dx 
^  f(x)  e  y 
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Thus   Sqn  41  becomes. 


d 
dx 


|f(x)dx 

e  y 


[f(x)dx 

=  g(x)e  , 


or 


r  f(x)dx 


|f(x)dx 

g(x)e  dx 


(42) 


For  our  case,    the    integrating  factor  is: 
dp^    Jl p . 


r  1- 


^3  M 


P 


g_ 


V/e    therefore   obtain: 
(Q  ^C,p)^2 

p 


^Vxv^V,) 


r 


P< 


P 


p. 


'  £3+1 

(CVC2P)^2 


PV| 


iy!^l  dp' 

dp'  I     ^ 


where 


p' =  a  so-called  "dummy  variable" 

p^=  air  density  at   beginning   of  altitude    interval. 
V/e    impose   the    condition  that   at    h^    (where    p   =po    ), 


V        =  V     =0.      Thus   v/e   have: 

XW  X 


£3 ,1 

(Ci^C2P)^2 
P 


M 


(V^x) 


p 


p' 


/P*i 


(Ci^CopT^ 
P' 


(3V2 


dVo/^ 


dp' 


dp.(V3) 
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Solving  for    (V^^  -  V^)    we   find; 


V  -V 
xw     X 


g 


^p 


p° 


^S-^-4 


LP' 


•  C4 


d  Vxw   ,    , 
7-    dp    . 


dp' 


where    C^    =   o     ? 


Vv'e    nov/  make    a    chan.cre   of    variables   as; 


^-^ 


Po 


(44) 


(45) 


The    new  limits   of   integration  become    1   and    ^     and  we   obtain: 


xw        X 


(Q-qp  ^ 


C3 
Co 


+1 


iQ^cp7 


7 


C4 

d  ?7  ^ 


(46) 


V/e   will    finally   be   interested    in  the    case  where 
V-,.,  has   a    linear  variation  v/ith  altitude.      For   the   present, 
we   shall    consiaer  the   more   general   case   where 

V^  ~-  a(h-hj  (47) 

\7e  must  express  V  as  a  function  of  ?7  for  Use  in 
Eqn  46.  V/riting  the  exponential  density  law  in  terms  of  the 
dummy  variable,      p'     ,    we    obtain: 

p'    _       -p(h'-h,) 


or 


Po 

h'-ho 


e 


9/p 


/3 


Thus,    we   have 


V, 


xw 


a 


In 


P>/p' 


n 


L  P 


(4S) 


14 


p 


which  may  be  rewritten  (since  "q"  ~  '?  ^  ^^^ 

In 


V      -  Cr 


P/p 


^  OC 


P    \ 

Differentiating,    v/e    find: 

n-1 


^  p/^ 


/^ 


(4ga) 


d  V 


xw 


^   in"  Wv     ■ 


(49) 


d?/  ft"  71 

Substituting  Eqn  U9   into   Eqn  46,    we   finally  obtain: 

C4 


/.^ 


Vxw-Vx 


Ci^C2P7 


^ 


"1         L^ 


iH 


(3"  7? 


n     p/t}  ^'^  ■        (50) 


This   equation  gives  an   integral  relation  for   V^^^  -  V^   for   any 
given   shoor   profile   of    the  form   shovm  in   Eqn  47. 

V/e    now   specialize   to    a    linear   profile   by   letting 
n    =  1   and   obtain: 


M<w-Vx 


a 
/3 


.^r  .^^1   1^^ 


Ci^CgP     1^2 


7 


di 


(51) 


which  may  be  rewritten  as 


OC 


/    1     -Cr,p 


n 


c. 


11 


(52) 


where: 


C 


C 


(k4,  )Vol 


Ci 


mc^  AP 


TR< 


Vol 


^3,1    __    1^K 


C2  "   R/R.-pK 


(j 


(53) 


C. 
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We  must  now  integrate  relation  (52).   Since  we  are 
considering  altitude  intervals  no  larger  than  500  ft,  the 
upper  limit  of  tlie  integr^.1  in  relation  (52)  is  quite  close 


to  unity.   For  example,  from  Ref  4  we  have:  ^3    =   a^  aaq 


and  therefore 


Pi 
P 


-  e 


h-ho 
23.000 


500 
23sOOO 


for  h 


h. 


500 


and   thus  we   have    Pr°  =   1.022.      Numerical   examini-tion  of  the 

y 

term  within  the   bracket   in  Eqn   52   reveals  that  this   quantity 
is   almost  a   linear   function  of  >;    in  the   range    >^    =  1  to 
'n    =  1.022.      Thus,    we   may  approximate  this   by  a   Taylor's 
series   about    7i     ^   1   with   all   but  the   first   tvro    terms   neglected, 
That    is,    we   may   express 


g(7) 


1 -C5P 


C3.1 

Co 


^ 


by: 


One  easily  finds  through  differentiation  that 


,'(>?  =  i)  ^  T-  \cl^^\ 


1  *C5p 


and  therefore 


g(  ^)  -1  " 


Mc.-v, 


Substituting   this   in  relation    (52)    we   find: 
(3  1      ^^2       h-CsP 


(T^-Dj      jll 


(54) 


1^: 


Since  7)  =   ^       over  the  range  of  integration  v;e 


may  replace  2J2  by  dv   and  obtain: 
V  -V  =  ^ 


^C4 


^'[^-S:*^lr?^p]«-^M  '^   '551 

This   may   be   integrated  directly  to   give: 


Vxw-Vx  = 


P^CQ+DP 

C2        1  *CspJ 

94"1 


i-[i-(§*^)i§^pk^) 


(56) 


or: 


^xw     -^ 


'OC 


X      /3(C,.l)[i.(|3.,)^4^pJ 


1  + 


[iig-'),4 


C5P 


2      n-CsPJlP 


p 


ro  _-] 


-1 


.(57) 


For  the   second   solution  we  assume  p    is   constant 
over  the  altitude   interval  and  write: 


^x      dh    dt  2dh 

Substituting  this  in  Eqn  30  and  rearranging  yields: 


(5^) 


dVx 

dh 


VxwVx 


V| 


C3P0-C1 
CsPo+Q 


(59) 


v/here      Po  =  sir  density  at   h  =  h      (lower  edge   of   shear  layer). 

H  \/ 

Subtracting       ,  ,^    from  both  sides   of  this  equation  and 
d  h 

rearranging  gives: 

d(VxO^Vx) 
dh 


C3P0-C1 


dV, 


9_/ w_w    )  -    J±_j 

Lc^po-qj  v.^^"^^""^'     cih 


xw 


(60) 
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Assximing  that   V     and  f>     are    constant,    we   find   the 
integrating  factor  to    be: 


Jc2po*C-l  "^    ^     ^    C2P0-C1    Vz2j        ^     ^C2p.*CiV,2 


(61) 


Thus    Eqn  60   becomes 
d 


dh 


C3P0-C1  g'^ 


=    e 


c.:tPo-ci  gh 

C2Po*Ci'^ 


d  V; 


xw 


dh 


(62) 


Once   again  assuming  a  wind   shear  given  by  V^^ 
or  (h   -   h    ),    v;e   obtain: 


o 


dh 


c^p.-q    gh 


C3P,-Ci     gh 


or,    after  integrating  , 

CtP«-Ci  gh 
e^^^i  ^  (V -Mc) 


h  =  h 


=  CX 


h=h<, 


C2P0*  C, 

[CaPo-q. 


2  C3p,-Ci  Jh 

Vz     _C2Po*Ci     Vj,^ 


h^h. 


h=ho 


at   h  =  hQ,    V^^  ~  "^x  "^    ^^xw  ~  ^x^o   ^^^  thus   we   can   solve   this 

to   obtain 

_    Oi  P"  -  Ci    g(h-ho) 


oc 


C3Po-qJ     g 


C3po-q    g(h-ho)  ^ 
1-  e~  C2p«*<^i        V22 


(63) 


which  may  be  written  as: 

_f(h.)(h-hj 
(Vx^V^)-(VxwVx)oe. 


f(ho)L 


1-  e 


-f(ho)(h-ho)' 


(64) 
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where 


f 


(hj-.C3£^^ 


C2P0.C1   v|     • 

Summarizing  these  two  solutions,  v/e  have: 
Solution  I:   Variable  density,  V  =  constant, 

a 


(65) 


VxwVx  =  p(c^.i)g(h) 


Solution  II:   Constant  density,  V  =  constant. 


1 


(66) 


(v^vJ-(\^Y<), 


-r(hJ(h-ho)  cc 


f(ho)  _ 


1  -  e 


-r(ho)(h-ho) 


(67) 


f  Ch  ^  -  ^3P°-ci  _9_ 
^''^^^  -  C2P0-C,    v| 


We   shall   now  show  that   Solution  II   is  approximately 
equal  to   Solution   I   for   small  altitude   intervals.      Expanding 
g(h)    and   using  the  definition  Cr   =  _^   ,    we   obtain: 

g(h)  =  2lS^Pzi%t2l£il' 
y  VC5P 

Therefore,    v;e   may  write: 


^  /    g  ,    C1-C3P 

p(Vi)g(h)  =  p(^a^i^-^T^    • 
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But      ^— ?  >>  1        and    GO    v/e   may  v;rite 


Thus    v/e   observe   that 


p(  C4^i)g(ho)  =  -f(ho) 


(68) 


Using  the  experimental  density  law  for  small 
altitude  intervals,  we  find: 


^°-1  -ft 


(h-ho)  , 


Thus,  v;e  have 


1.g(ho)(#-1)  =  1  ■^g(ho)(3(h-ho)  . 


Also,  for  h  -  h^  small,  this  represents  the  first  two  terms 
for  the  Taylor  expansion  of  e 
Thus  we  may  write: 


g(ho)p(h-ho) 
_^g(ho)(h-ho) 


Or,  for  small  altitude  intervals  we  have: 


Po 


Using   relation    (68)    we   have 
"l^q(ho)(-^°-1  ) 


C4+I  7^g(h,)(h-h,)(C4^1) 

=    e 


^0  ^•lC4^1 


4  ■  '  -f(ho)(h-ho) 

=     e 


Thus,    Solution   I   becomes: 


xw  > 


(X 

^       f(ho) 


1  -  e 


-f  (ho)(h-ho) 
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which  is  identical  to  Solution  II  if  (V,   -  V,  )   =0. 

Thus,  we  see  that  Solution  II  is  of  sufficient 
accuracy  if  h  -  h^  is  sufficiently  small,  i.e.,  if 
(h  -  h^)  ^  500  ft. 

C«        Comparison  of  Analytical  Solution  with  Computer 
Solution 

As  a  check  on  the  accuracy  of  Eqn  6?  (Solution  II), 
a  comparison  was  made  with  computer  solutions  of  Ref  2. 
These  comparisons  were  made  at  large  and  small  altitudes  as 
well  as  large  and  small  wind  gradients.   Figures  1,  2,  and  3 
present  the  solutions  from  Ref  2  as  compared  with  Eqn  6?. 
The  solid  curve  was  obtained  from  Ref  2  and  represents  a 
computer  solution.   The  circled  points  v;ere  obtained  from 
Eqn  67  vdth  the  apparent  mass  neglected  (K  =0).   The  dashed 
curve  was  obtained  from  Eqn  6?  with  the  apparent  mass  factor 
(K)  assumed  to  be  0.5. 

Comparison  of  the  circled  points  with  the  solid 
curve  indicates  that  Eqn  6?  is  extremely  accurate  for  a  11 
altitudes  and  wind  gradients.   It  is  seen  also  that  the 
v;ind  response  error  is  significantly  increased  due  to  the 
apparent  mass  characteristics  of  the  balloon,  consequently 
indicating  that  the  a  pparent  mass  m.ust  not  be  neglected  in 
problems  of  this  nature. 
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FIG  1.  WIND  SHEAR  AND  BALLOON  VELOCITY 
(WITH  APPARENT  MASS  FACTOR  OF  0.0 
AND  0.5)  AS  A  FUNCTION  OF  ALTITUDE 
(h^=  15,000  ft,  ex  :=  0.14) 
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ho  =  15,000  ft.  (4,570  meters) 

a  =  WIND   GRADIENT  =  0,28 

DIA.  =  6  ft.  (1.83  meters  ) 

WT.  =  330  gms. 

GAS  -  HELIUM 

V7  =  13.9  meters /sec.  (45.6  ft /sec.) 
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FIG  2.  WIND  SHEAR  AND  BALLOON  VELOCITY 

(WITH  APPARENT  MASS  FACTOR  OF  0.0 
AND  0.5)  AS  A  FUNCTION  OF  ALTITUDE 
(h^=15,000ft.,   a  =  0.28) 
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FIG  3.  WIND  SHEAR  AND  BALLOON   VELOCITY 
(WITH  APPARENT  MASS  FACTOR  OF  0.0 
AND  0.5)  AS  A  FUNCTION  OF  ALTITUDE 
(h^=55p00  ft,   oc:=0.28) 
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III.  i>;':7sU:Iu:m^nt  of  the  ai-p,::-j:t 

KhSS   OF  I:   SPHFI^F-LIKE  BALLOCI^ 

A  ,         Ex po rime nta 1  Pro c ed n  ro 

Iv'hen  a  body  moving  in  a  fluid  experiencos  an 
accelerating  force,  the  x'caction  to  this  force  is  such 
that  the  chan£^,e  in  velocity  with  time  is  slov/er  than  expected, 
appearing  as  if  the  body  had  a  mass  greater  than  the  actual 
measured  value.   This  experience  is  explained  as  the  effect 
of  the  so-called  "apparent  n-ii-ss",  v;hich  has  been  calculated 
from  potential  flow  equai,ions  for  several  relatively  simple 
bodies.   For  several  more  complicated  parachute-like  bodies, 
the  apparent  mass  has  then  been  established  experimentally 
(Ref  5). 

A  spherical  balloon,  used  to  measure  v/inds  as  it 
rises  through  the  atmosphere,  experiences  accelerations  and 
decelerations  as  it  passes  through  the  layers  of  horizontal 
winds.   The  trajectory  calculations  of  the  preceeding  chapter 
show  that  the  apparent  mass  is  of  significant  influence  and 
m.ust  be  considered  in  any  determination  of  the  v/ind  pattern 
by  m.eans  of  an  ascending  balloon.   Therefore,  the  objective 
of  the  following  study  is  to  establish  experimentally  the 
apparent  m.ass  for  a  wind-detecting  balloon  fabricated  from 
flat  gores  and  having  an  approximately  spherical  shape  v;hen 
properly  inflated. 
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To   ir.cr.suro   t,  lie   cpparcnt    m.:-~G    cf  -the-    bulloonr., 
experimental    pro  ccdurc;T.  v;ere  adopted  :.  r,  used    in   r.^rcchute 
tec!!nolof';y    (Ref    5).      i.odiryin^^    t:ie    referenced    rr.etnod    for  this 
purpoe-Cj    r-he    balloon.':   v/erc   drop^^cd    from  a    hei^iit    of   about 
25    fo>-t    in    a    h£.n<;^ar-likc    labor,  tory   at    the    Unlvercity   of 
••"linnesott  .      A    knovm  v.'oi^.ht   v/as    sucpended   beneath  the   balloon 
and   allov;ed    to    strike   the    floor.      This   impact    provided   an 
instantineous   d'eceleration   of   the    balloon   system  v/hile    the 
dra^   of   the    system  momentarily   rem.ained    the    seme    as   before 
the    impact.      From  this    test    procedure,    the   apparent   mass 
can   be   determined   in   the  lollov/in^   manner. 

The   forces   actinr;;   on  the    balloon  during   the   fall 
are   shov;n   in  Fig  4.      If   the   balloon   is    in   steady   state 
(dv/ct   =  0)    at   this   tim.e,    Nev-'^on's   second    la\-/  provides  the 
follo-.;ing   relation: 


or 


W^  +  Wg  +  W^   +  Wg-  B  -  D    =   0 

D    =  W^  ^  W(3  +Wl  ^  Wg  -  B     • 

V/hen  the    lost   v;eight  V/^   strikes   the  floor,    the 


ecuation   of  m.otion   is 


W,  -W^  -W„  -  B  -  D   ::  [Wa^Wg^Wb    ^ 

A  o  D  n 


m 


dv 
dt 


Solving  for  ra'  and  considering  the  deceleration,  dv/dt,  to 


(69) 
'(70) 


(71) 
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DRAG  (D) 


WEIGHT  OF  ACCELEROMETER  (Wa) 


LOST  WEIGHT    (  W  ) 


FLOOR 


FIG  4.    FORCES    ACTING    ON    THF   EXPERIMENTAL 
SYSTEM    DURING    DESCENT 


;-,-.  ---^.-ativc,  v;g  ootam 


oc 


m 


(72) 


by  r  to  ciir.plif'y  Icte-'  calculc.tio 


nz .    one  obtains; 
(73) 


In  the  preceding  analysis  it  v;as  assur...d  that  the 
balloon  did  not  accelerate  during  its  fall,   Hov;ovcr,  in  the 
actual  case  if  was  found  that  the  balloon  v;as  still  accel- 
erated v.'hen  the  lost  weight  hit  the  floor.   The  eo_uations  of 
motion  were  modified  accordingly.   If  we  let  a^  "be  the  accel- 
eration of  the  balloon  just  previous  to  the  impact  of  the 
lost  weight,  the  equation  of  motion  becomes: 


Wa  +W^+W^+Wg-B-D  = 


WA^\A(5^WL^WB^rn' 

g 


(74) 


The  net  v/eight  of  the  balloon,  F,  is  given  by: 

This  is  a  balance  reading  when  the  balloon  system  is  set 
upon  it,  and  hence  is  one  of  the  measurements  made  during 
the  testy   If  Eqn  75  is  substituted  into  hqn  74,  v;e  obtain 


F-D  = 


'"  Wa  ^Wg  ^Ws^Wl  ^  ^' 


ao 


(76) 


Using  the    nev;  definition  and   after  the   v;eight    strikes 
the   floor,    vie  have 

r  "1 

d  V 


F  -  D  -W,     = 


Wa  ^Wg  ^Wp       ^  ^/ 

T       III 


dt 


n--.' 
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Solvinr-  for  ;:i'  ,  v/o  obcair. 


F  -  D  -  Wl   Wa  ^  Wr.  ^  vy 


dv/dt  g  ^^  ^ 

Or,  multiplyinr  by  ii^ 

The  right-hand  side  of  Eqn  79  is  easily  deter- 
mined, v.'ich  the^  exception  of  the  dra^-^,  D,   Velocity  raeas- 
uremants  v.^ere  made  i-r.c  ohc  drag  could  bo  calculated  from 

D  ^^pV^^CoS   .  (SO) 

Hov.'evir,  it  v;as  found  that  the  Reynolds  number  of  the 

experiricnts  v/as  close  to  the  tr:  nsitional  ro.nz^   of  a  sphere, 

(3.45  ::  TO-"  and  1.77  x  10^  for  the  2.0  meter  bailee,  i;  2.60  x  10^ 

and  "i  .31  x  10  for  the  2.44  meter  balloon),  and  the  0,-^   values 

became  soraev;hat  uncertain. 

In  order  to  avoid  detailed  experiments  to  evaluate 

precise  Cn's  for  these  models  at  the  operating  Reynolds  numibers, 

a  reiterative  aooroximation  orocess  was  used.   Since  a   is   

...  o     

alv;c.ys  small  (0.004  ^  ^q  "^  0.022  g's)  one  m^ay  assum:e  that 
the  balloon  is  in  steady  state  before  im.pact  hence  drajr  is  khDv;n 
and  './e  may  calculate  v;'.   Usinr  this  value  anc"  IJev/ton's  Second 
Lav.'.  ::e   obtain  from  Sen  7h   a  now  dra^  term; 

y 

Using  tnis  value  for  the  drag  and  substituting  it  into  2qn  79? 
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one   obc^i;-i;   c    ncv;  v:  luo    .'or  the   C-pycr-^nt   v;ei,"ht„      Continued 
reitoru-cions   betwee-n  ^'.qciz   79   ^"i^-    ■'■"i    v;ere    v)crrorr.ioa   until 
the   vi-lue   of  '.7'    ch?.r.;;;ed    in  e£.ch   ;.r-ocesc   by   less   than  0.5/&. 


3.      ::c.';oi 


The  experimental  syste:-.  for  apparent  maos  measure- 
ments consisted  of  a  bsllocn,  balloon  rigging  and  v/eights, 
acceleroraeter,  and  a  velocity  nac:  suring  device  (Fig  5). 

1  )  'Balloons 

Tv;o  helium-filled  balloons,  constructea  of  2  mil, 
mylar,  v/c  'o  used  in  the  -^escs.   Tl-c  lorger  ballooa  v;as  2.44 
meters  in  oiameter  (8.0  feet)  and  the  smaller  balloon  had  a 
diaaeter  of  2.0  meters  (6, $6  feet).   The  2.44  meter  balloon 
(Fig  6)  had  24  gores  and  vras  nearly  spherical.   The  2.0  meter 
balloon  (Fig  7)  had  only  12  gores  and  obvious  deviations 
fro."  spherical  shape  v;ere  observed  (  See  Fig  S).   The  volume 
of  anis  balloon  v;as  caj  culated  assuming  the  balloon  was 
composed  of  several  special  ungula s  of  a  right  circular  cyl- 
inder. 

In  prelim.inary  vrark,  the  balloons  were  filled  to 
a  specific  inflation  pressure  -.;i::h  air.   With  this  arrange- 
m.entj  the  m.ass  of  the  G;'Gte!r:  becir.ie  extremely  large  and  since 
the  mass  of  the  included  gas  must  be  considered  in  the  calc- 
ulations of  acceleration,  it  appeared  necessary  to  reduce 
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FIG    7     PHOTOGRAPH       OF 
METER       BALLOON 


INFLATED  2 

(12    GORES) 
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FIG  8.  INFLATED    12  GORE     BALLOON    SHOWING 
NONSPHERICITY    OF   CROSS -SECTION 
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the  system  mass  by  introducing  a  lighter  gas  into  the  balloon. 

Using  helium,  the  balloons  developed  a  free  lift 
of  $i  lbs  for  the  2.0  meter  balloon  and  9^  lbs  for  the  larger 
one.   Then  using  a  lost  weight  of  six(6)  and  tentlOO  lbs, 
respectively,  the  resultant  dovmv/ard  force,  the  system's  mass, 
and  its  dre.g  were  considerably  decreased.   This  greatly 
simplified  test  procedure,  calculations,  and  reduced  the 
probability  of  acceleroraeter  damage  due  to  excessive  loads 
and  rebound  acceleration. 

Using  this  system,  the  lost  weirht  could  easily 
be  varified,  causing  impact  velocities  and  Reynolds  number  to 
vary  in  the  desired  manner. 

2)   Ballon  Rigging  and  V^eights 

The  balloon  rigging  consisted  of  a  harness  of 
three  0.017  inch  steel  cables,  equispaced  over  the  balloon 
surface,  fastened  together  at  the  top  and  meeting  at  a 
confluence  point  below  the  balloon  and  supporting  the  "lost" 
weight  (Fig  5). 

3)   Accelerometer 

A  Statham  accelerometer  with  a  range  of  *  3  g's 
was  used  to  measure  all  accelerations  of  the  balloon  system. 
The  accelerometer  was  fastened  to  an  aluminum  sphere  section 
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of  14  inch  chord  length  which  was  attached  to  the  base  of 
the  balloon.   This  arrangement  formed  a  rigid  system  of 
balloon,  plate,  and  accelerometer ,   In  preliminary  work  the 
accelerometer  was  suspended  by  a  single  cord  a  few  inches 
below  the  balloon.   This  suspension  was  unsatisfactory 
because  the  accelerometer  became  excited  from  vibrations 
of  the  balloon  when  the  lost  weight  contacted  the  floor  and 
the  static  and  dynamic  equilibrium  of  the  balloon  was  disturbed, 

4)   Velocity  Indication 

The  velocity  of  the  descending  balloon  had  to 
be  measured  prior  to  ground  impact.   For  this  purpose  the 
following  method  was  employed:   A  copper  screen  was  stretched 
^  inch  above  an  aluminum  plate  and  both  surfaces  were  placed 
on  the  floor.   The  plate  and  the  screen  were  connected  to 
a  Century  oscillograph  which  recorded  the  instant  of  contact 
between  the  screen  and  the  aluminum  plate.   The  time  base 
of  the  Century  read-out  system  was  used  to  measure  the  in- 
terval between  the  instants  at  which  a  small  steel  ball,  sus- 
pended beneath  the  main  part  of  the  lost  weight,  and  the 
main  portion  of  the  lost  weight  itself  contacted  the  wire 
screen. 

The  ball  was  heavy  enough  to  depress  the  screen 
upon  impact,  but  light  enough  to  allow  the  screen  to  break 
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contact   v/ith  the    ground    plate.      The    lost    weight    provided 
the    permanent   deflection  shown    in  the   oscillogram. 


4)      Data   Recordi 


ng 


The  accelerations  and  velocity  indications  were 
recorded  and  a  sample  trace  is  shown  in  Fig  9.   Note  that 
the  balloon  acceleration, is  averaged  over  a  finite  time 
period.   This  is  permissible  since  the  aerodynamic  drag, 
which  changes  with  the  velocity,  is  small  compared  with  the 
other  terms  of  the  system  shown  in  Eqn  69. 

C.        Results 

A  summary  of  results  is  presented  in  Table  1  and 
Fig  10.   The  value  of  K,  defined  as 

is  predicted  by  potential  flow  theory  to  be  O.5OO  fof-  a 
sphere  (Ref  3).   As  is  shown  in  the  table,  the  experimental 
values  are  quite  close  to  this  analytically  derived  value. 
Reviewing  the  results,  one-notices;  in  the  table  as  well,  as  in  the 
graph  a  trend 'Of  Increasing ■ apparent -mass  with  decreasing  Rey- 
nolds number  .   This  may  be  explained  as  being  a  consequence 
of  the  varying  flow  pattern.   It  was  mentioned  that  the 
2.0  meter  balloon  deviates  more  from  the  spherical  form 
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FIG  9   SAMPLE    OSCILLOGRAPH    TRACE 


SPHERE 
DIAMETER   (ft) 

6.56            6.56 
(2   r-ETER) 

8                8 
(2.44  JAETER) 

SPHERE 
VOLUME    (ft.) 

ia.2 

U1.2 

268.1 

268.1 

REYNOLDS    NO. 

345,000 

177,000 

260,000 

131,000 

IMPACT 

VELOCITY  (ft/sec.) 

7.1 

4.5 

5.4 

2.7 

AVERAGE  IMPACT 
DECELERATION 
(q's) 

0.54fi 

0.503 

0.467 

■0.461 

AVERAGE   DESCENT 

ACCELERATION 
(q's) 

0.022 

0.004 

0.022 

C    ,  DRAG 

COEFFICIENT^ 

0.119 

0.182 

0.074 

APPARENT  MASS 
BUOYANCY 

0.519 

0.550 

0.485 

0.512 

X    ESTIMATED 


TABLE  I.  EXPERIMENTAL  RESULTS  AT  VARIOUS 
REYNOLDS  NUMBERS  FOR  A  FREELY 
DESCENDING    SPHERE 
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than  the  larger  one.   This  difference  in  the  shape  could 
account  for  the  7%   difference  in  the  apparent  mass  factor 
K  for  the  two  balloons.   Since  K  for  a  cube  is  O.67O 
(Refs  5  and  6), and  the  projection  of  the  2.0  meter  model 
is  more  polygonal,  an  app^-ii^ent  nass  factor  higher  than  the 
one  for  a  sphere  appears  to  be  reasonable. 
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IV.   CONSLUSIONS 

The  trajectory  caXcul^-tions  in  Section  II  illus- 
trate the  response  error  and  the  pertinent  equations  show 
the  terns  v;hich  essentially  cause  the  lack  of  motion  of  the 
balloon  behind  the  actual  ;vind.   Results  shov;  that  the 
apparent  mass  is  a  significant  factor  and  must  be  included 
in  the  analysis  of  such  a  system  (Figs  1  through  3).   For 
this  analysis,  a  theoretically  derived  value  of  K  =  0.5  v/as 
used  as  the  apparent  mass  factor. 

The  trajectory  c£;  Iculations  described  in  Section  II 
compare  favorably  with  computer  solutions.   The  more  exact 
calculations  as  well  as  the  solution  of  the  equation  of  motion 
appear  to  be  acceptable  for  use  in  data  evaluation. 

An  analysis  of  the  wind  response  equations  (Eqn  No 
6?)  shows  that  a  minimum  error  can  be  obtained  if  the 
following  criteria  are  considered: 

1)  The  b£lloon  should  rise  as  slov/ly  as  possible. 

2)  The  apparent  mass  factor  K  is  of  importance 
and  must  be  kept  as  small  as  possible. 

Also  examination  of  relation  (14)  indicates  that; 

3)  The  ratio  of  the  weight  to  drag,  W/Cy,S, 
should  be  a  minimum. 

In  addition  to  these  considerations  are  several 
other  items  such  as  size  for  radar  soundings  and  balloon 
strength  and  shape.   The  pursuit  of  these  circumstances 
exceeds  the  scope  of  this  study. 

The  apparent  mass  factor  used  in  the  calculations 
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was  assumed  to  be  K  =  0.5.   Experiments  shov/ed  that  the 
apparent  mass  factor  of  actual  balloons  is  close  to  the 
theoretical  value.   A  certain  dependency  (Table  I,  Fig  10) 
of  the  apparent  mass  factor  upon  Reynolds  number  is  indicated, 
but  under  the  scope  of  the  present  program  this  phenomenon 
was  not  investigated  further. 
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